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Abstract 

We consider a wide class of determinants whose entries are moments of the so-called semiclassical func- 
tionals and we show that they are tau functions for an appropriate isomonodromic family which depends 
on the parameters of the symbols for the functionals. This shows that the vanishing of the tau- function 
for those systems is the obstruction to the solvability of a Riemann-Hilbert problem associated to certain 
classes of (multiple) orthogonal polynomials. The determinants include Hankel, Toplitz and shifted- 
Toplitz determinants as well as determinants of bimoment functionals and the determinants arising in 
the study of multiple orthogonality. Some of these determinants appear also as partition functions of 
random matrix models, including an instance of a two-matrix model. 
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1 Introduction 

The connection between orthogonal polynomials and monodromy preserving deformations of differential 
equations has been observed repeatedly [25, 18, 19, 16, 20, 5, 6]. It can be used to produce special solution 
to Painleve equations by choosing appropriately the measure of orthogonality and its dependence on 
deformation parameters [25]. In the eighties the notion of isomonodromic tau function was developed 
[23, 21, 22] for systems of monodromy preserving deformations of rational connections (under some 
genericity assumptions); the philosophical attitude was that -in the context of Birkhoff-Riemann-Hilbert 
problems- the tau-function should play the same role as the classical Theta function plays for the Jacobi 
inversion problem on algebraic curves. Indeed the vanishing of a Theta- function on the Abel map of a 
divisor determines its "specialty", namely the (im) possibility of finding suitable meromorphic functions 
with prescribed divisor. Similarly the vanishing of the isomonodromic tau function indicates that a 
Riemann-Hilbert problem is not solvable. This was proved first for the Schlesinger isomonodromic 
equations by Malgrange [26, 27, 28] and more recently for connections with irregular singularities by 
Palmer [33] (although under some genericity assumptions). 

In this perspective a connection between Hankel determinants and isomonodromic tau functions is 
not surprising since the vanishing of the Hankel determinant of the moments of a measure determines 
whether the orthogonal polynomials exist or not and -on the other hand- their existence is equivalent to 
the solvability of a RHP [19]. 

The precise functional relationship between these (and similar) determinants and isomonodromic tau 
functions was established in [4, 5] for the Hankel case and in [20, 6] for shifted Toeplitz determinants of 
so-called semiclassical moment junctionals [25, 29, 30]. 

In a seemingly distant area, moment determinants appear as partition functions for random matrix 
models [31, 32], which are known to be tau functions (in a different sense) for the KP, Toda or 2Toda 
hierarchies, depending on the case. For example the partition function of the Hermitean matrix model 
defined as 



where M is a Hermitean matrix and dM is the Lebesgue measure on the vector space Ti n of Hermitean 
matrices of size n x n; here V{x) represents any scalar function, but for the sake of the example we can 
take it to be a polynomial. 
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By means of Dyson-Mehta theorem [32] it is known that Z^ is (up to an inessential normalization 
constant) the same as a Hankel determinant for the moments of the measure (or weight) e~ v ^ dx 

Z x n MM oc A n := detf^cK^xn-i , Hj := j aPe~ v ^ dx (1.2) 
In addition the same quantity can be expressed as a multiple integral over the spectra 

Z\ MM oc J dx x ...J dx n A 2 (X)e-^U v (^ , A(X):=H(x j -x k ) (1.3) 

j<k 

This correspondence between matrix models and Hankel determinants is much deeper, inasmuch as 
the correlation functions for the random eigenvalues can all be described in terms of the orthogonal 
polynomials for the above measure. 

The relationship between Z\ MM and isomonodromic tau functions was established in [5, 4] but in 
this paper we will provide a new proof which applies to more general cases as well. An example is the 
following two-matrix model, recently introduced in [8] ; it consists of pairs M\ , Mi of positive-definite 
Hermitean matrices of size n x n with a partition function defined as follows 

Z 2MM / / dM dM ! ( L4 ) 

J J det(Mi+M 2 )« V ! 

It can be shown that this partition function is also a moment determinant [8] 

?2MM A , rTl - f f e-W- y ^> dxdy 3 k 

Z n := det[X] ifc , 2 jk := / / ■ x J y (1.5) 

JM.+ jr + x + y 

Similarly to the Hermitean-matrix model case, also here the spectral statistics can be solved in terms of 
certain biorthogonal polynomials for the pairing induced by the integrand above [7, 8] 

A third class of moment determinants that we will consider here and were considered in [6] are the 
shifted Toeplitz determinants 

A f n := det[/^ + j-_fc]o<j,fc<n-l (1-6) 

for moments defined as in (1.2). These enter as tau-functions for various integrable lattices of Toda/Ablowitz- 
Ladik type [15, 14, 13]. In the ordinary Toeplitz case [t = 0) they too can be interpreted as partition 
functions for unitary matrix models. 

A last class of multi-weight moment determinant we will consider here arise in the theory of multiple 
Pade approximation for Markov functions [1] . Let e Ve ^ dx , j = 1, . . . , K be a collection of (semi classical) 
weights with polynomial V^'s (for the sake of this example). For a point n = [ni, . . . ,njc] in the lattice 



we define the following determinant of size n = \n\ := Ylt=i 
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u {1) 
r«i-l 



u {1) 



Pn+m -2 



Ast := det 



~7 WY 
M 



u {K) 



rn-1 



(JC) 
L i n+n K ~2 



^ := xh~ V ^dx 



(1.7) 



Remark 1.1 VFe ignore whether there is a matrix model related to this determinant. 

The point of the mater is that for all these determinants, the dependence on the coefficients of the 
potentials V is isomonodromic (in the sense of Jimbo-Miwa-Ueno) and that they can all be identified 
with the isomonodromic tau function of [23]. The paper is devoted to this relationship. 

1.1 Organization of the paper and outline of the proof 

We first review very quickly the setup of [23, 21, 22] in Sec. 2, in order also to establish the notation 
and some results that are used later. We will need to compute explicitly the isomonodromic tau function 
for a very simple (almost trivial) class of isomonodromic families of upper-triangular matrices (Sec. 2.1) 
since this result will be used repeatedly. In Sec. 3 we recall in more detail the notion of semiclassical 
moment functional and associated Markov (Weyl) function. 

The last four sections are devoted to establishing the precise relationship between determinants of 
type Hankel (Sec. 4), shifted- Toeplitz (Sec. 5), Cauchy/bimoment (Sec. 6) and multiple-weight (Sec. 7). 

We briefly comment on the simple idea of proof since it is the same in all cases with changes only in 
the details of the implementation. 

The first step is to consider the relevant (multi)orthogonal polynomials associated to the determinant 
under inspection and characterize them by means of a Riemann-Hilbert problem. Most part of this step 
is actually well known (for example [18, 19] for the case of orthogonal polynomials) and can be found 
-for the different cases- in the literature. 

The second step is to establish an ODE whose (generalized) monodromy is constant by construction 
under deformations of the weights; this is accomplished by transforming the RHP into an equivalent one 
with constant jumps. This step requires the weights to be semiclassical in order to obtain a rational ODE 
and thus to fall within the framework of isomonodromic deformations. At this point the isomonodromic 
tau function recalled in Sec. 2 is constructed. We then recognize that the changes in the degrees (or multi- 
degrees) of the polynomials are suitable Schlesinger transformations and thus we can prove that the ratio 
of the corresponding tau function is the same as the ratio of the corresponding moment determinants. 

Therefore to complete the description of the relationship it is sufficient to identify r and the deter- 
minant for a "trivial" case (the "initial condition" ) . This is always the case of polynomials of degree 
and determinants of size (which equal 1 by convention). The isomonodromic family corresponding to 
degree-zero (multi) orthogonal polynomials is -in all cases- an instance of the "trivial" class of Sec. 2.1 
for which the isomonodromic tau function can be computed easily and explicitly. 
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2 Isomonodromic deformations: a (short) review 



The content of this section is lifted (with adapted notation) from [23, 21, 22] and we refer the reader 
ibidem for more details. 

Consider an rank-r ODE in the complex domain 

V(x) = D(x)W(x) , D{x) £ Mat rxr (2.1) 

with an r x r matrix D(x) whose coefficients are rational functions. We will denote by A the divisor of 
singular points (including oo if the case) and by a G A a point in it; for a £ A we denote by ( = £. a the 
local parameter at a, namely £ a = x — a if a ^ oo or = - for a = oo. We will denote by d a + 1 the 
order of the pole of D{x) dx at x = a. 

If the differential D{x) dx has a pole at a of order > 2 then the point is called an irregular singu- 
larity 3 , otherwise it is called a Fuchsian singularity. 

Suppose that near a singularity x = a we have the following (possibly sectorial) asymptotic expansion 

* Or) = G(l + *iC + >2C 2 + • • -)e T(C) (2.2) 
where T(£) is a diagonal matrix of the form 

da 

T(C) = y £T ] Ci + T ln( (2.3) 
j'=i 

It is understood that T = T (a) , y = Y. , , G = G( a ) are matrices that depend on the particular singularity 
under exam; when necessary we will use a subscript to remind of this fact. We recall the following fact: 

Proposition 2.1 ([35]) A sufficient condition for the expansion (2.2) to be valid in suitable sectors 
around an irregular singularity x = a is that the leading coefficient of the singularity of D{x) dx is 
ad-regular namely it has simple spectrum (i.e. distinct eigenvalues) . 

We stress -which is crucial for our setting- that the above condition is only sufficient and not necessary; 
in fact some of our cases below will violate the condition but still retain the validity of the expansion 
(2.2). 

The expansion (2.2) is in general only asymptotic, namely the series V(C) is only formal; a theorem 
[35] shows however that for each (sufficiently small) sector there is a bona fide solution whose asymptotic 
expansion matches (2.2). Two solutions ^f, having the same asymptotics in different sectors are related 
by a right constant multiplier 'I'(x) = ^f(x)S, called a Stokes matrix. Given two singularities a, a £ A 
and suitable sectors near them, the right multiplier for the solutions Vl/, Vl/ with the given asymptotics is 
called connection matrix. 

Additionally, analytic continuation of any solution to the punctured plane C \ A induces a (anti) 
representation of the fundamental group called the monodromy representation for the connection 
(d x -D(x))dx. 

We will assume that oo is a singularity and choose the normalization of Goo = 1 (this can always be 
accomplished by a left-multiplication of by constant invertible matrix). 

3 When counting poles for D(x) dx we must recall that at x = oo the differential dx has a second-order pole. 
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Definition 2.1 The generalized monodromy data consist of the totality of the Stokes and connection 
matrices, the monodromy representation and the Tq^'s appearing in (2.2). 

Remark 2.1 This definition is redundant since the data as defined above are not entirely independent 
or "minimal" . This is of no consequence since we are going to study only transformations that preserve 
those data. 

Definition 2.2 The isomonodromic times are the coefficients Tj A , j > 1 Va G A as well as the 
location of the singularities. 

In [23] it was shown how to deform the isomonodromic times while preserving the generalized monodromy 
data of the connection. 

Definition 2.3 The isomonodromic tau function r is a function of the isomonodromic times (de- 
pending parametrically on the generalized monodromy) defined via the following first-order PDE 



5lnr := £ resl* (y^Y^x)^*)] dx (2.4) 



where the symbol 5 stands for the total differential in all isomonodromic times. 

Of course closedness of the differential (2.4) needs to be proved, and that was accomplished in [23]. 

We point out that in [23] they were working on the assumption that all singularities of D{x) met 
the sufficient condition of Prop. 2.1; however it can be seen -going through the proofs- that the only 
information that they used was the validity of the expansion (2.2) with a diagonal T^{Q) near each 
singularity a £ A 

2.1 "Trivial" isomonodromic systems 

We will need the explicit expression for r in a very simple situation. Of course the simplest situation 
would correspond to a diagonal connection, for the problem would be immediately solvable. Suppose 
that D{x) is upper-triangular and also that at all singularities a £ A we have 

*(x)~Y (a) (x)e T («) (x) (2.5) 

where T,Jx) is diagonal and Y,Jx) upper triangular and (formally) analytic in the local parameter 
C = C.. Denoting by Yd{x) the diagonal part of Y the isomonodromic tau function is simply given by 

dlnr = resTr(y- 1 yjdT) (2.6) 



For definiteness let 

V{{x) * 

D(x) 



Vi(x) 



(2.7) 



6 



with Vj{x) some rational functions. Then the solution has the following form 



-Vi 



(2i 



Near a singularity x = a we thus have that T(x) = — diag(Vi, . . . , V r ) s i ng is the singular part of the 
potentials at x = a whereas Yd(x) = exp ^— diag [Vj, . . . , Vr\ re ^j where the subscript denotes the regular 
part. 

Thus we have 



51nr = ^ res[V-)a )r eg5(Vj) a[ 



sing 



(2.9) 



where the symbol res denotes the sum of the residues at all singularities. 
A 

A simple exercise shows that 



lnr = ^res(V j )a,re a (V-)a, 



3=1 



sing ■ 



(2.10) 



It should be noted that there could be other isomonodromic parameters in the upper triangular terms 
of D(x) (and hence of Y(x)), but clearly the tau function does not depend on them. 



2.2 Schlesinger transformations 

The family of monodromy-preserving transformations can be slightly enlarged by allowing discrete trans- 
formations that change the entries of To, (a) by integers 



-o,(a) 



! 0,(a) 



+ L 



(a) 



L 



(a) 



diag(^ a)i i...^ a)ir 



-(«).} 



(2.11) 



subject to the constraint YlaeA ^^V) = 0- These transformations can be achieved by multiplying on 
the left the solution \& by a rational matrix R{x). Such matrix can be algorithmically constructed from 
the entries of the matrices {C (a) , K ,} ae ^ as explained in [21]. They are called Schlesinger transfor- 
mations and they do not change (by definition) neither the monodromy representation nor the Stokes 
and connection matrices. Schlesinger transformations form a discrete Abelian group (a lattice) of trans- 
formations generated by the so-called elementary Schlesinger transformations. These are such that 
the matrices {L (a) } ae _4 have only two nonzero entries, a +1 and a —1. For example one may have 
L ia) = diagl 1 * 0, . . .) and L (oo) = diag(0, 0, -1,0, . . .). 



The details can be found in [21], but to give a taste of the type of computation we consider the 
elementary Schlesinger transformations where £(00) = diag(0, . . . , 0, 1, 0, . . . , 0, —1,0, . . .) = <7jj, where 
the 1 is on the i-th position and the —1 is on the j-th. 

Denote Y( 00 )(x) = Y(x) = 1 + Y^jL\^i x ~ J '> the problem is then to find a second formal series 
Y{z) = 1 + X]j=i YjX 3 and a (polynomial) matrix R(x) such that 



Y(x)x a » = R{x)Y(x) 



(2.12) 
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The solution is straightforward but tedious: R(x) must be linear in x and of the form R{x) 
Matching the coefficients of the powers of z carefully one can find R$ and R\ 



Rq + xR\. 



R(x) 



1 

1 


a\ 

Q>i—1 








hi... bi-i 


x + A 


b\ ... bj-i 


h 


bj+i ...b r 




a-i+i 
a 7 _i 


1 

1 








1 

6,- 











a j+1 
a r 






1 

1 



(2.13) 



where 



A :-- 



Y~, 



-Y, 



(2.14) 



It appears quite obviously that 

Proposition 2.2 The problem above admits solution if and only ifY± d . ^ 0. 
The importance for us lies in the following theorem -which we paraphrase- 



Theorem 2.1 (Thm. 4.1 in [21]) Given two connections D{x)dx and D(x)dx related by an elemen- 
tary Schlesinger transformation with a +1 in the i-th diagonal element of and a — 1 in the j-th one 
of Lr a i\ then the corresponding isomonodromic tau functions r, r are related by 



5 In 



where 



H 



( Y (a),l)ij 
{G(a'))ij 

(G (4 G w)^ 



SH 



if a = a' 

if a = oo, a' ^ oo 
if a ^ co, a' = oo 

if a, a' / oo, a' ^ a 



(2.15) 



(2.16) 



With these preparations we are ready to investigate the relationship between several types of moment 
determinants and isomonodromic tau functions. 
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3 Markov functions for semiclassical moment functionals 



We will consider a rather wide class of determinants of matrices whose entries are moments of (collections 
of) weights. 

The weights we are considering are all of the semiclassical type as defined in [25, 29, 30, 5, 11]. This 
means that they are of the form /j,(x) = e~ v ^ with V'(x) an arbitrary rational function. 

They are integrated over contours jj which can be arbitrary contours in the complex plane as long 
as all integrals x k [i(x) dx are convergent integrals. The range for k will be either N or Z, depending 
on the situation; a detailed description of the contours can be found in [5, 6] and we refer thereto for a 
more detailed discussion. 

We will choose arbitrary complex constants x~ for each contour jj and use the notation 

/ x k fi(x) dx := Xj / x k fi(x) dx = fik (3.1) 

We will also use the notation n : C — > C to indicate the locally constant function that takes the (constant) 
value Kj on the corresponding contour 7,-. The Markov function (sometimes referred to as Weyl function) 
for these semiclassical weights is simply defined as the locally analytic function on C \ U7j given by 

W{x) := [ , k := £ x iX7j (x) (3.2) 

J " ^ X j 

At times we will denote by k (by abuse of notation) the support of k. The function W{x) has logarithmic 
growth at the hard— edges, namely endpoints of contours jj where [x is 0(1). In this case one verifies 
that W{x) = 0(ln \x — a\), where a is the hard-edge point. 



4 Hankel determinants: orthogonal polynomials 

4.1 Riemann-Hilbert problem of orthogonal polynomial 

Let {p n (x)\n = 0, 1, 2, ...} be the (monic) OPs that satisfy the following orthogonality 

Pn{x)p m {x)e~ v{x ' ) dx = h n S nm . 



(4.1) 



They can be uniquely characterized by the Riemann-Hilbert Problem (RHP) described hereafter. 
Define for z G C \ M+ the matrix 



T(z) :=T n {z) :-- 



p n (z) C\p n ](z) 



C[p]{z) 



1 

2i7T 



p(x)e v ( x ' dx 



x — z 



(4.2) 



The above matrix has the following jump-relations and asymptotic behavior that uniquely characterize 
it [18, 19, 16, 17] (we drop the explicit dependence on n for brevity) 



r + (x) = r_(x) 



1 xje- v ^ 



, x € Yj , 



r(*)~ i + E 



z n 
z~ n 



(4.3) 
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Near a hard-edge z = a (where -we recall- V(z) is analytic) we have T(z) = [0(1), 0(ln \z — a\)]: more 
precisely we have 

1 



T(z) = G a {z) 



G a (z) = 0(1) 



(4.4) 



Replacing the orthogonality condition (4.1) by the above jump (and boundary) conditions (4.3, 4.4) we 
obtain the Riemann-Hilbert problem for the OPs. 



4.2 A problem with constant jumps: ODE 

If we introduce 

^ := re"^ 2 ^ 3 
then we the RHP (4.1) is turned into an equivalent RHP 

^ ~ (1 + 0(z- 1 ))z na3 e-^ a3 



1 Kj 

1 



(4.5) 



(4.6) 



It is understood that near all singularities of V(z) the matrix ^(z) has the singularity implied by the 
multiplication by exp(— \gz) an d near the hard-edges it has the singularities implied by (4.4). 

Proposition 4.1 The matrix ^f n (x) solves a rational ODE ^f' n ( x ) = D n (x)^f n (x) where D n {x) is a 
matrix with rational coefficients with the same singularities as V'(x) and with simple poles at the hard- 
edges with nilpotent residue. 

Proof. Noticing that det ^ = det Y = 1 we have that ^'(z)^~ 1 (z) has no jumps across the contours Tj 
and hence we have an ODE 

9'(z) = D n (z)9(z) (4.7) 

where D n {z) is analytic in the plane punctured at the singularities of V and at the hard-edges. 

Remark 4.1 The matrix D n {z) can be written explicitly in terms of the recurrence relations as in [5]. 

Near a singularity z = a of V' we have ^ n { z ) = G a {z)e~^' 73 with G a (z) a formally invertible asymptotic 
series in z — a (or 1/z if c = oo, in which case = G a (z)z nas e~ 2" CT3 ) and thus 

V\z) i 



D n {z) 



G(zy 1 a 3 G(z) + 0(l) 



Near a hard-edge we have 



* n (z) = G a (z) 



W Hz - a) 
1 



1 



and thus 



D n (z) = G~ L (a) 



2m(z~a) 





G a = 0(l) 



G a (a) + 0(1) 



(4.8) 



(4.9) 



(4.10) 



Using Liouville's theorem it is seen that D{z) (which depends on n) has the same singularity structure 
of V'(z) and additionally some simple poles with nilpotent residue at the so-called hard-edges (if any). 
Q.E.D. 
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Definition 4.1 We define 



A., 



det[/jL i+ j]o<i,j< n -i, (J>j ■= / x j e v ^ dx 



(4.11) 



It is known that A n is -up to a ra-dependent constant- the partition function of the Hermitean matrix 
model [31]. 

Theorem 4.1 The Jimbo-Miwa-Ueno isomonodromic tau function r n associated to the isomonodromic 
family ^ n can be normalized so that 

T n = A n r (4.12) 



with tq given by (2.10) for V\ 



-Vo 



■\V. 



Proof. We show that T n+ i/r n = A n+ i/A n . 

We first observe that \P n +i and ^ n are related by an elementary Schlesinger transformation as in 
Sect. 2.2, namely 



ra+l 



X C n Cl n 

± 



(4.13) 



where c n and a n are uniquely determined in terms of the matrices Yg that appear in (4.3). The relation 
(4.13) is nothing but a rephrasing of the three-term recurrence relation 



Xp n (x) = p n +l(x) + C n p n {x) + a n p n -i{x) . 

From the representation (4.2) we have that 

1 



(4.14) 



Y\ 12 = coefficient of z n 1 in 



1 f p n (x)e v ^ dx 
2iK 



x — z 2in 
Since the monic polynomials admit the determinantal representation 



z~ n ~ X / p n {x)x n e- v(x) dx (4.15) 



Pn(x) = -r— det 



Mo Mi 
Mi ■•• 



Mn-1 • • • 
1 X 



Mn 
Mn+1 



M2n-1 



(4.16) 



we see immediately that f^p n (x)x n e v ^ dx = A ^ +1 . 

By Thm. 2.1 and formula (2.16) we can choose the normalizations of r n+ \ and r n so that 



T n +1 



-2ot(Ki) 



A n 



12 



+ 1 



A, 



(4.17) 



This immediately implies that r n = ToA n (since -by definition- Aq = 1) 
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The isomonodromic function tq is associated to the isomonodromic family It is straightforward 
to verify that the solution of the RHP (4.3) for n = is given by 



r (*) 



2in J >c x—z 

1 



*o(*) := 



_Y_W e 2 f e V{x, dx 

e 2 / 

2m x — z 

e 2 



(4.18) 



A straightforward computation yields for Dq(x) 



V'(z) 







V\z) 



--: D (x) 



(4.19) 



where the star denotes an expression which is irrelevant to the present considerations (but can be easily 
computed). We are thus in the setting of Sec. 2.1, and thus To is given by formula (2.10). Q.E.D. 

5 Shifted Toeplitz determinants: biorthogonal Laurent and Szego poly- 
nomials. 



We assume that the moments jij are defined for all j £ Z and introduce the notation 



A< := det 



/ H He+i • • • Ht+n-i \ 
Ht-\ Hi " ' He+n-2 

He~n+l He~n+2 ■ ■ ■ He 



I He He+i 
He-i He 



T n {x) := det 



He-n+i He-n+2 

V 1 X 



He+n \ 
He+n-i 

He+i 



(5.1) 



Ag = l, Ai =0 



It is apparent from the determinantal expression that is "orthogonal" to all powers x 1 »•• 

Define the monic polynomials as 

x n + 0(x n ~ l ) 



e ._ i J n _ 

Pn ■ a e 

L - X n 



Consider the matrix 



We have the RHP 



C[f] 



1 [ f(z)n(z)dz 



2m 



z — X 



j 2*7r(-l) 



n-l 



A* 



(5.2) 

(5.3) 
(5.4) 



r, =r_ 



1 K/J,(z) 

1 



(1 + Oix- 1 ) 



x n 
x~ £ - 



r ~ < 



(5.5) 
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where 



Note that 



Lt,, 



(~)"A 
W 

A e 



t+i 



-A* 

n+l 



A £ 



detGl = n+1 n ~,\^ n n = 1 



e \2 



det T e n (x) 



(5.6) 
(5.7) 



5.1 Constant-jump RHP and ODE 

We now postulate ^(x) = e~ v ^ with V'(x) a rational function. We introduce the new matrix 



1 K 

1 



(5.8) 



* = (i + o^- 1 )) 



x 

x 



-i-i 



e 2 



* ={l + 0{x)) 



1 

x r 







e 2 



This matrix solves a similar RHP with constant jumps and hence also an ODE 

d x ¥ n = Dl{x)¥ n 

As in Prop. 4.1 one can see that the connection has the same singularities as V'(x), where 

V(x) := V(x) - (n - £ - 1) Inx . 



(5.9) 



(5.10) 



(5.11) 



plus nilpotent-residue simple poles at the hard-edges. 

The various shifts n — > n ± 1 and I i— > t ± 1 all correspond to Schlesinger transformations (of more 
general form than the one used for OP). In [6] we considered (from a different approach) the shifts 



(n,l) 
(n,£) 
(n,t) 



(n + l,£) 
(n + l,£+l) 
(n,£+l) 



Circle move 
Line move 
Circle-to-line move 



(5.12) 



Of course the "Line move" is simply the composition of the other two; the reason for the naming was 
that a circle move implies that the OP satisfy a recurrence of the form 



(5.13) 



for some constants indicated anonymously by the *. This recurrence is typical of the Szego orthogonal 
polynomials on the circle, whence the name. 

Likewise a line move corresponds to a recurrence of the form 



Xpi = P^+l + +Pi + 



(5.14) 



typical of the recurrence relation for orthogonal polynomials on the real line. The circle-to-line move is 
yet a different recurrence that intertwines the two. 
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Remark 5.1 These mixed recurrence relations were important in the study of integrable lattices [15, 13] 
and put on the same grounds the integration of all known lattices (Toda, relativistic Toda, Ablowitz-Ladik, 
Volterra) and also some unnamed generalizations. 

Remark 5.2 If all the moves are Circle-moves, then the resulting OPs are Szego polynomials on the unit 
circle (if the weight is a real weight on the circle, that is, or generalizations thereof). If all the moves are 
Line-moves then the polynomials are the usual (non Hermitean) orthogonal polynomials. 

Following Thm. 2.1 and formula 2.16 together with the expressions (5.6) one can verify that -in all cases- 

A e ' T 1 ' 

Indeed a Schlesinger transformation that changes the exponents only at infinity (like in Sect. 4) corre- 
sponds to a Line- Move and r'/r is given by the (1,2) entry of the Y\ term in the expansion at infinity. 
The other moves are Schlesinger transforms that change one exponent at oo and another exponent of 
formal monodromy at x = 0. In particular (referring to (5.6)) 

1. for a Circle Move, ^f- oc (G e n ) 12 ) 

£ 

2. for a Circle-to-Line Move, -ppx oc (G„)ii. 

Once more oc means up to any expression that does not depend on the isomonodromic parameters. 
Inspection of G„ yields (5.15). 

The simplest way to completely describe the relation between and is to note that the shifted 
Toplitz determinant is (up to a reshuffling of rows) a Hankel determinant of the same size for the weight 
x n ~^fi(x) dx. This corresponds to considering the isomonodromic equation for 

*i(x) := <(x)x^ (5.16) 

¥ n (x) + = %(z). 



1 K 

1 

+ Oix- 1 )) z^e"^ * ¥ n = G e n (1 + 0{x)) (x)ff3 (5.17) 



which shows that -as the reader may have observed- the polynomials <p l n and Q^-x are si m ply the 
(unnormalized) orthogonal polynomials for the weight x n ~ i ~ 1 fi(x) dx (on the same contours k). 
This implies that 

r l n = r Al (5.18) 

where to is the same expression in (2.10) but with V\ = — V% = — \V . This result reproduces the results 
of [6] for shifted Topliz determinants and [20] for (unshifted) Toplitz determinants with special symbol. 
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6 Bimoment determinants: Cauchy biorthogonal polynomials 



Let a(x),(3(y) be two arbitrary semiclassical weights and consider the bi- moment functional 

,a(x)f3(y) dx dy 



I - •- 



x % y ] - 



x + y 



(6.1) 



where the symbol k XiV denotes any linear combination of contours as long as the integrals make sense. It 
is understood that the contours of integration for the weights a and (5 must be such that the intersections 
of the contours k x H (— K y ) consist of isolated points, so that the integral is well defined. 
We denote by A n the principal minors of size n 



A n :— det[/]o<jj<„_i 



(6.2) 



The biorthogonal polynomials are two sequences of (monic) polynomials {p n (x)}n£N, {ln(y)}n£N 
of exact degree n such that 



a{x)(3{y)dxdy _ ZW 



x + y 



A, 



(6.3) 



(Note that the form of the constant appearing in the RHS is actually a small theorem). 
They admit the determinantal representation 



Pn(x) = ~r— det 



00 



'-On 



I-n—10 ■ • ■ In— In 



, Qn(y) 



A, 



det 



2", 



oo 



1q n-l 



1 



T T i i r - 

-LnO • • • -Lnn—l y 



(6.4) 



Remark 6.1 If a, (3 are positive measures and the contour of integration are the positive real axis then 
it was shown in [8] that the bimoment matrix is totally positive; moreover the zeroes of the OPs are 
real, simple and interlaced. It was also shown ibidem the these BOPs satisfy a four-term recurrence 
relation^ . 

In [7] the BOPs and some auxiliary polynomials were introduced in order to describe the Christoffel- 
Darboux identities, which in turn have important applications to the study of the spectral statistics for 
the Cauchy matrix model [8]. In the present work we will only need the following characterization for 
the sequence of (monic) OPs p n (x) (with an entirely specular characterization for the q n (y))- 

Proposition 6.1 (Prop. 8.2 in [7]) Consider the Riemann-Hilbert problem (RHP n ) of finding a ma- 
trix T(z) = T n {z) such that 

1. T{z) is analytic on C \ (k x U (— K y )) 

4 It was proven for positive weights on the real axis, but the proof is entirely algebraic and hence applies without changes 
to this semiclassical situation. 
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2. T(z) satisfies the jump conditions 

?(*)+ = ?(*)- 

r(*)+ = r(*)- 



1 K x a(z) 

1 

1 

1 

1 Ky(3*(z) 







5. its asymptotic behavior at z = oo 9(2) / is 



1 



Z S — K. 



(6.5) 
(6.6) 



/?*(z) :=0{-z) 



T(z) =1 + 



z n 
1 
0^ 



(6.7) 



Then such a T(z) uniquely characterizes the polynomial p n ( z ) fl s its (11) -entry. Moreover T(z) can 
equivalently be written as: 



T(z) 



PO,n Pl,n P2,n 
P0,n-1 Pl,n-1 P2,n-1 
P0,n-1 Pl,n-1 P2,n-1 



(6.8) 



where the first column consists of polynomials of the indicated degree and, for any function f(z) we have 
denoted 



h(z) := f(z) 



fi(z) := 



1 



f(x)a(x) dx 



1 



2iir L x — z 2iir L z + y 



h(-y)f3(y)dy 



(6.9) 



Remark 6.2 The detailed description for the entries ofT (although for a different normalization) can 
be found in [7]. 

Since in [7] the proposition was stated only for positive measures on the real axis, we briefly sketch the 
proof that the above RHP characterizes the OP's p n . The jump- relations and the asymptotics at infinity 
imply that the first column is an entire function bounded by z n and hence -by Liouville's theorem- a 
polynomial. Moreover the asymptotic at 00 imply that the (11) entry is a monic polynomial of degree n 
and the remaining entries some polynomials of lesser degree. 

The jump relations imply that the second and third column are obtained from the first by the 
expressions (6.9) (using Sokhotsky-Plemelj formula). Finally, the (33) entry of the asymptotic at infinity 
implies that 

j a(x)/3(y)dxdy 

p n (x)y J =0, j<n-l (6.10) 

x + y 

which is the orthogonality requirement. 
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6.1 Constant-jump RHP and ODE 



A simple transformation brings the above RHP into one with constant jumps. We will write a{x) 
e - v i( x ) an( j f3*(x) := (3(—x) = e~ V2 ^ x \ We introduce the new matrix 



V(z) :=T(z) 



exp i 



exp (^) 



exp 



(6.11) 



which solves a RHP with constant jumps 



*(z) 



*(z)- 



1 K x 

1 
1 



1 
1/% 
1 



(6.12) 



In the same way as for the OP we can prove 

Proposition 6.2 The matrix = ^f n ( z ) solves an ODE of the form 

*' n (z) = D n {z)^ n {z) 



(6.13) 



where the 3x3 matrix D n {z) has the same singularities as V{(z), V^{z) plus simple poles (with nilpotent 
residues) at the hard-edges. 



Proposition 6.3 The ratio of two consecutive tau functions satisfy 

Tn+l _ 



(6.14) 



Proof. By inspection of (6.7) we observe that the shift Vl/ n \— > ^ n +i is a Schlesinger transformation as in 
Sect. 2.2 involving the exponents of formal monodromy at infinity in the entries (1, 1) and (3, 3). Once 
again, using Thm. 2.1 and formula (2.16) we have (with a convenient normalization of tau-functions) 



Tn+l 



(6.15) 



It is immediate to see that (^1)13 is the coefficient in front of z n 1 in the asymptotic expansion of P2,; 
namely 



P2,n( Z ) 



a{x)P{y)dxdy _ x 

Pn\X) , . 1 . — { ) z 

(z + y)(y + x) 



Pn(x)y' 



,a(x)(3(y)dx dy 



(-)»* 



y + x 



+ 0(z 



-n-2^ 



=(6.16) 



+ 0{z- n -\%\7) 



This concludes the proof. Q.E.D. 
We thus prove 
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Theorem 6.1 We can choose the normalization for the Jimbo-Miwa-Ueno tau functions r n for the 
isomonodromic family ^ n so that it is related to the bimoment determinant A n (6.2) by 



Tn = ToA n 



(6.18) 



with r given by (2.10) with V x = V 2 = Y ^ , V 3 = (here the V h2 ,3 on the LHS refer 

to the Vj 's appearing in (2.10), while V\ 2 ° n the RHS refer to the potentials used in the definition of the 
weights). 

Proof. By Prop. 6.3 we already know r n = ToA„, so it only remains to compute To- The RHPo for 
is easily solvable and yields 



1 



1 f a dx 1 

X — Z 47T 5 



f f 



a(x)f3(y) dx dy 
(z+y)(y+x) 



1 C /3dy 
2in 



r ei 

Jk v y- 
1 



exp 



2V 1 +V 2 \ 
3 I 



exp(^) 



exp 



Thus the ODE d z — Dq{z) has an upper-triangular form 

_2Vi+Vi 

D (z 



YtYi 



3 



(6.19) 



(6.20) 



where the precise expressions for the upper-triangular entries is presently irrelevant. We are thus in the 
setting of Sec. 2.1 and we can use eq. 2.10. Q.E.D. 

Remark 6.3 We remark -in passing - that the spectral curve of Dq (i.e. its characteristic polynomial) 
is a completely factorized rational curve 



(6.21) 



7 Multi-measure moment matrices: multiple orthogonal polynomials 

Multiple orthogonal polynomials appear in the theory of simultaneous Pade approximations (for a review 
and interesting applications to number-theory and spectral theory of banded matrices see [1]). 

We will give a (unfairly) quick definition which allows us to get to the point as directly as possible. 
Consider K measures ^\x). Denote by n E a multi-component vector and |n| = Ylf—i n j [34]. 

The multiple orthogonal polynomials of type II are polynomials of degree \n\ defined by 



P H (x)x i n U \x) dx = , W : < £ < nj - l,j = 1, 



K. 



(7.1) 
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The system is called perfect if P^ has exact degree n := \n\. We can then assume that P^ is monic and 
the following determinantal expression are immediately deduced 



Pn(x) = -r- dot 



U W 



u {l) 
u {1) 

H'n+l 



u {1) 



IF 



u {K) 



u {K) 



(K) 



1 X 



X 



Ar 



det 



u (1) 



u (1) 



Mn+m -2 



^0 



u [K) 



(A") 
L l n+n K -2 



(7.2) 



They are characterized as the (1,1) entry of the solution of the following Riemann-Hilbert problem 



r n(^)+ = r n (x). 



1 

1 



(10 



--:T n {x)_M{x) 



T n {x) = {l + 0{x- 1 )) 



x 



-hi 



X 



-n K 



(7.3) 



We observe that the first row of the solution T = is written explicitly as 



Tn = P n (x) , Tij 



1 



(7.4) 



7.1 Constant jump RHP and ODE 

As it stands the RHP (7.3) cannot be related directly to the theory of isomonodromic deformations, since 
the latter refers to deformations of rational ODEs. For general weights the RHP cannot be reduced 
to a rational ODE. 

In order to reduce the problem to one with constant jumps it is sufficient to find a matrix G{x) such 
that G(x)Z 1 M (x)G(x) + is (locally) constant. 

The properties of G(x) depend on the class of weights we are considering. We will only consider 
semiclassical weights but there are situations (which we will consider elsewhere) in which the weights are 
not semiclassical and yet a more refined approach can be used to obtain a rational ODE. 
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7.1.1 Semiclassical weights 

If the weights are semiclassical and of the form fi£ = e~ Ve we can chose G(x) as follows 



G(x) = diag 
K 



(7.5) 



Vo := Yl V o > % " Vj = Vk ■ (7.6) 

3=1 fc^3 

We thus define 

* = tf* = T n (x)G(x) (7.7) 

which now solves a constant-jump RHP and following the same steps as Prop. 4.1 we obtain an ODE of 
the form 

%(x) = D H (x)V n (x) (7.8) 

where -D^(x) is a rational matrix with the same singularity structure as the (totality of the) Vjs as well 
as nilpotent-residue simple poles at the hard-edges (if any) . 

Theorem 7.1 The multi-measure determinant /S.^ in eq. (7.2) and the Jimbo-Miwa-Ueno tau function 
Tft for the isomonodromic family ^ are related by 

Tft = A h t$ , (7.9) 

with given by formula (2.10) with V\ = —Vq and Ve = Vt , £ = 2, . . . K + 1. 

Proof. The matrix ^ has the asymptotics 

tt*(x) = (l + — + ^ + ..\ x T -G{x) 
\ x x z J 

T H = diag(n, -n%, . . . , -n K ) . (7.10) 

Suppose n' = ft + a, where eg is an standard basis vector. The transformation ft i— > n' correspond 
to a Schlesinger transformation with (in the notation of eq. (2.12)) Gx\+i- Then Thm. 2.1 and formula 
(2.16) state that 



5 = y wn ( 7 - n ) 

i<n. 



and inspection of the matrix yields 



Y la4+l <x [ P H (x)x n ?^Xx)dxK^- (7.12) 

where the proportionality is within signs or other constant factors that are irrelevant because of the 
inherent ambiguity in the definition of isomonodromic tau function. 
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With this we have proved that the ratio any two adjacent multi-measure moment determinants is the 
same as the ratio of the corresponding isomonodromic tau functions 



and thus 



A- ~~ T- 

L - l n 1 n 



/ — i n 1 o ' n 



Once more, the computation of is simple because is an upper-triangular connection 

1 f //%)d£ 



1 Wi(a:) ...Wk(a;) 

1 



G{x) , := 



★ * ★ 
V{(x) 



V K (x) J 

Once more we are in the setting of Sec. 2.1 and formula (2.10) applies. Q.E.D. 



(7.13) 
(7.14) 

(7.15) 
(7.16) 



8 Conclusion 

It is apparent that the proofs for the four classes of moment determinants are essentially a repetition of the 
same idea. In order not to overburden the paper we did not consider other cases that nevertheless might 
have interesting applications. These include (but are not necessarily limited to) the cases of Angelescu 
and Nikishin systems for semiclassical weights. In fact the case of Sec. 6 is an instance of Nikishin system 
[1] but we singled it out because of its independent appearance in the context of random matrix theory 

[8]- . . \ 

Another important multimatrix model that relates to biorthogonal polynomials is the Itzykson-Zuber 

matrix model with a partition function 

Z™ := J dMi dM 2 e- Wl ( Ml '- W2 ( M2 ' +TrMlM2 (8.1) 
which corresponds to the following bimoment determinant [12, 2] 

cx detl jk , l jk = f f dxdyx j y k e- Vlix) ~ V2(y)+xy . (8.2) 



This determinant can also be related to multiple orthogonal polynomials [24] but for non semi-classical 
weights. While the general idea of the proof Z^ z is an isomonodromic tau function is the same, the 
details are substantially different. The main reason is that the relevant isomonodromic family does not 
have an asymptotic expansion as in (2.2) at oo and thus the definition of isomonodromic tau function 
of the Japanese school cannot be applied as is. In particular the leading coefficient of the relevant ODE 
has a very degenerate spectrum and the asymptotics involves fractional powers of z [3, 10]. The proof in 
this case will be detailed in a forthcoming paper [9]. 
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